OU-TAP-36 
May 1996 



SCALE-INVARIANT CORRELATION FUNCTIONS OF 
COSMOLOGICAL DENSITY FLUCTUATIONS IN THE 
STRONG CLUSTERING REGIME 

Tailiei Yano and Naoteru Gouda 
Department of Earth and Space Science, Graduate School of Science, Osaka University 

Toyonaka, Osaka 560, Japan 
El-maih yano, gouda@vega.ess.sci.osaka-u.ac.jp 



ABSTRACT 



We have investigated the scale-invariant solutions of the BBGKY equations 
for spatial correlation functions of cosmological density fluctuations and the 
mean relative peculiar velocity in the strongly nonlinear regime. 

It is found that the solutions for the mean relative physical velocity depend 
on the three-point spatial correlation function and the skewness of the velocity 
fields. We find that the stable condition in which the mean relative physical 
velocity vanishes on the virialized regions is satisfied only under the assumptions 
which Davis & Peebles took in his paper. It is found, however, that their 
assumptions may not be general in real. 

The power index of the two-point correlation function in the strongly 
nonlinear regime depends on the mean relative peculiar velocity, the three-point 
correlation function and the skewness. If the self-similar solutions exist, the 
power index in the strongly nonlinear regime is related to the power index of the 
initial power spectrum and its relation depends on the three-point correlation 
function and the skewness through the mean relative peculiar velocity. 

Furthermore it is found that the mean relative physical velocity should 
have the values between and the Hubble expanding velocity. When the mean 
relative physical velocity equals to the Hubble expanding velocity, there might 
exist self-similar solutions in which the power index of the two-point correlation 
function in the strongly nonlinear regime is independent of the initial power 
index n. 



Subject headings: cosmology: theory-large scale structures-correlation function 



- 2 - 



1. INTRODUCTION 

In our universe, there are many kinds of structures such as galaxies, clusters of galaxies 
and super-clusters of galaxies. It is not known clearly still today how galaxies and the 
large-scale structures have been formed after the universe was born. This is one of the most 
important problems in the cosmology. In the standard scenario of the large-scale structure 
formations, it is considered that small fluctuations at the early time grow as time increases 
due to the gravitational instability. Hence it is very important to clarify the evolution of 
density fluctuations. Here we consider the density fluctuations of the coUisionless particles 
such as dark matters because our interest is concentrated on the effect of the self-gravity. 
When the density fluctuations are much smaller than unity at the early time, time evolution 
of the small density fluctuations can be analyzed by making use of the linear theory. In 
this regime, we can understand analytically how the small fluctuations grow(Peebles 1980, 
1993). When the amphtude of these small fluctuations becomes as large as unity, that 
is, the fluctuations become to stay at the quasi-linear stage, we can not make use of the 
linear theory. In this regime, the higher order perturbation methods or the Zel'dovich 
approximation is often used for the analysis. Moreover the density fluctuations continue to 
grow as time increases and then the amplitude of these fluctuations becomes much larger 
than unity. At last the caustics of density fields appear anywhere in this regime. In this 
strongly nonlinear regime, the analytical approach is impossible. However the nonlinear 
phenomena of the self-gravity is very interesting and mportant for not only the large-scale 
structure formations, but also the academic interest on the nonlinear dynamics. Hence we 
believe that it is very necessary to understand clearly the nonlinear behavior of the density 
fluctuations. For example, we are interested in the two-point spatial correlation functions in 
the strongly nonlinear regime. In this regime, it is found from the A'"-body simulations that 
the two-point correlation functions obey the power law. This result is reasonable because 
the self-gravity is scale-free. Then the power index of the two-point correlation function is a 
good indicator representing the nonlinear dynamics of the self-gravity in this regime. And 
the power index is related to the clustering pattern of coUisionless matters. Hence it is very 
important to study what physical processes determine the power-index of the two-point 
correlation in the strongly nonlinear regime. Moreover it is interesting to analyze whether 
the power index in this regime depends on the initial conditions or not. 

These problems have been usually analyzed by iV-body simulations(Frenk, White, & 
Davis 1983; Davis et.al. 1985; Suto 1993 and references therein). The method by iV-body 
simulations is straightforward for tracking the evolution of the density fluctuations in the 
nonlinear regime. In the simulations, however, high spatial resolution is necessary for 
estimating the correlation functions and the mean relative peculiar velocity on very small 
scales in the strongly nonlinear regime with good accuracy(Jain 1995). But at the present 
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time, the computer ability has not enough high resolution. Then we believe that the 
dynamics in the strongly nonlinear regime have never been completely verified by using the 
A^-body simulations until now. 

There are other methods for the analysis of the nonlinear density fiuctuations. One 
of them is the analysis by the BBGKY equations. The work by Davis & Peebles (1977; 
hereafter DP) is a pioneer for the analysis by the BBGKY equations. They showed the 
existence of the self-similar solutions for correlation functions under some assumptions. 
Then it is shown that the power index 7 of the two-point spatial correlation function in 
the strongly nonlinear regime is related to the initial power index n of the initial power 
spectrum P{k) as follows: 



This result is very interesting, but it also seems strange because in general systems 
might forget the initial memories in the strong nonlinear regime due to the nonlinearity of 
the self-gravity. One of the assumptions that DP adopted is called the stable condition. 
This condition is that the mean relative physical velocity in the strongly nonlinear regime 
is equal to zero. This condition was tested by A'"-body simulations (Efstathiou et al. 1988; 
Jain 1995) and this condition is not completely verified. Furthermore the stability of the 
self-invariant solution in the strongly nonlinear regime, which DP derived, was investigated 
by the perturbation theory (Ruamsuwan & Fry 1992) and it is found that the solution is 
marginally stable. 

As for the physical process determining the power index in the strongly nonlinear 
regime, there are other analysis besides one proposed by DP. One of them is given by 
Saslaw(1980). He concluded that the power index 7 approaches to 2 by using the cosmic 
energy equation under some assumptions while some numerical simulations do not support 
this result (Frenk, White & Davis 1983; Davis et al. 1985; Fry & Melott 1985). 

There is another idea as follows; when the initial power spectrum has the sharp cut-off 
or the initial power spectrum is scale-free with negative and small initial power index, then 
there appear anywhere caustics of the density fields. In these cases, the power index is 
irrespective of the detailed initial conditions after the first appearance of caustics on the 
small scales around the typical size of the thickness of caustics(in two-dimensional systems, 
they correspond to the filament structures of highly clustered matters) . The power index is 
determined by the type of the singularity which is classified in accord with the catastrophe 
theory. This idea is verified in the one-dimensional system(Kotok & Shandarin 1988; 
Gouda & Nakamura 1988, 1989), the spherically symmetric systems(Gouda 1989) and the 
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two-dimensional systems (Gouda 1996). In these cases, it is suggested that 7 f=i on the 
small scales. 

As we can see from the above arguments, we believe that there are uncertainties about 
the physical processes which determine the value of the power index. In this paper, we 
examine the conditions which determine the power index by analysing the scale-invariant 
solutions of the cosmic BBGKY equations. The analysis of the nonlinear clustering in 
the strongly nonlinear regime by the BBGKY equations has advantage rather than that 
by the //-body simulations. This is because the BBGKY equations directly deal with the 
statistical quantities such as correlation functions and this analysis is free from the artificial 
coUisionality due to finite numbers of particles which might appear on the small scales in 
the A^-body simulations. However there is technically difficulty in solving the BBGKY 
equations; in our analysis, the BBGKY equations are translated to the moment equations 
by integrating the equations over velocity for convenience. Then these moment equations 
for time evolution of the A^-point spatial correlation function includes the terms including 
A^ + 1-point spatial correlation functions. Furthermore A^-th moment equation includes 
the terms including A^ + 1-th moment. In general, these equations have infinite hierarchy 
and can not be closed at the lower order spatial correlation and moment. Hence we should 
close the equations in taking some assumptions. For example, DP took the following 
assumptions; the three-point spatial correlation function C is represented by the products 
of the two-point spatial correlation functions ^ as follows: 



where Q is a constant. Some observations suggest that this relation with Q ~ 1 holds at 
^ ~ 1. Furthermore DP assumed that the skewness of the velocity fields vanishes. Adding 
one more assumption, DP closed the BBGKY equations. As mentioned before, DP used the 
stable condition in deriving the power index 7 given by eq.(l). In this paper, we reexamine 
the scale-invariant solutions of the BBGKY equations and estimate the value of the power 
index and its dependence of the initial power index n when the above assumptions and 
stable condition are changed. And we analyze whether the stable condition is satisfied 
or not and how the mean relative peculiar velocity would be in real. Indeed, Jain(1995) 
claimed that the stable condition has not yet been verified by A"-body simulations. 

Furthermore we investigate whether there is possibility that the power index in 
the strongly nonlinear regime does not depend on the initial power index n even if 
the self-similarity of the solutions is satisfied. Recently Padmanabhan(1995) suggested 
the possibility that the power index is independent of n by using the pair conservation 
equations. 
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In this paper, we examine the above problems by using the cosmological BBGKY 
equation in the strongly nonlinear regime. 

In §2, the cosmic BBGKY equations are briefly reviewed according to DP and 
Ruamsuwan & Pry (1992). We show the scale-invariant (power-law) solutions in the strongly 
nonlinear regime and their properties in §3. Pinally, we devote §4 to conclusions and 
discussions. 



2. BASIC EQUATIONS 

In this section, we briefly review the derivation of the cosmological BBGKY equations 
according to DP and Ruamsuwan & Pry (1992). 



2.1. Cosmological BBGKY Equations 

Here we derive the BBGKY equations from the ensemble mean of the Vlasov equation 
in the expanding homogeneous and isotropic background universe. In this paper, we 
consider the only Einstein-de Sitter universe because we are interested in the scale-invariant 
solutions of the correlation functions and the self-similarity of the solutions and so it is 
necessary that the background universe has scale-free. The A^-body correlation function 
is the statistical quantity which is given by the ensemble mean of the iV-products of 
the one-body distribution functions. Then the BBGKY equations can be derived by the 
ensemble mean of the Vlasov equation. The Vlasov equation for the one-body distribution 
function f{x,p) is given by 

dt mo?dx"' ^dx'^dp"' ' ^ ' 



0(a;i) = / -d-'x2(rp2, (4) 

a J \X2 — a^il 



where m is the mass of a particle, a is the scale factor and G is the gravitational constant. 
The ensemble mean of the Vlasov equation is 

W^J^M_w^^)^0. (5) 
dt ma? dxi dxf dp" 

We define the statistical functions which are given by the ensemble mean as follows: 



6(1) = {f{x^,Pi)), 



(6) 
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P2(l,2) = (/(a;i,pi)/(a;2,P2)), (7) 
P3(l,2,3) = (/(a;i,pi)/(a;2,P2)/(a;3,P3)), (8) 



where 6(1) is a function of only the momentum because of the homogeneity in the 
background universe. The ensemble mean of the Vlasov equation (5) is rewritten by using 
above functions(eqs.[6] ~ [7]) as follows: 



at a opi J \x2 — xi\'^ 



This is the first BBGKY equation. As we can see from eq.(9), the time evolution of the 
one-body distribution function depends on the two-body correlation function. 

The following is the second BBGKY equation for the two-body correlation function: 

5P2(1,2) d 



dt dt 



{f{Xi,Pi)f(x2,P2)) 



dp2{l,2) Gw? d 



f^p3(l,2,3)ci^X3d^P3+(l-2), (10) 
ma^ oxi a opi J x^i 

(11) 



-''31 _ -^3 -^1 
3^31 ~ 1^3 - Xi\^' 



We can get the /^-body correlation function by the same way as follows 

9p,3(1.2.3) d 



dt dt 
ap4(l,2,3,4) ^ d_ 

dt dt 



{f{Xi,Pi)f{X2,P2)f{X3,Ps)), (12) 
{fiXi,Pi)f{x2,P2)f{X3, P3)f{X4, Pa)) , (13) 



In our analysis, however, we use the only second BBGKY equation as shown later. 

The second BBGKY equation is the time evolution equation of the two-body correlation 
function. In this equation, the three-body correlation function is involved. In general, the 
time evolution of the A'^-body correlation function depends on the N + 1-body correlation 
function. 

We define the following irreducible correlation functions, c and d: 

P2(l,2) = 6(1)6(2) + c(l, 2), (14) 
P3(l,2,3) = 6(l)6(2)6(3) + 6(l)c(2,3)+6(2)c(3,l)+6(3)c(l,2) + d(l,2,3). (15) 
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Here b, c and d mean 

- b{pi), (16) 

c{i,j) = c{xi,pi,Xj,pj), (17) 

d{ij,k) = d{xi,Pi,Xj,Pj,Xk,Pk)- (18) 

The first and the second BBGKY equations are rewritten by using the above functions, 
respectively, as follows: 

db{l) Grri^ d r x^^ 



dt 



+ / ^c(l, 2)d^X2d% = (first BBGKY) (19) 

a opi J X21 



dc(1.2) p'f dc(l,2) Gm^a^l) f x^, o^ ,s 

+ 2 a + ^c{2,3)d^xsd'ps 

at ma^ oxf a opf J X31 

/ ^'^(l' 2, ?,)d^x^d^p^ + (1 ^ 2) = (second BBGKY) (20) 



+ 



2.2. Velocity Moment 

We are interested in the power index of the two-point spatial correlation function in 
the strongly non-linear regime. Hence the equation which we use in our analysis is the 
second BBGKY equation. Moreover we use the velocity moment equations which are given 
by multiplying the second BBGKY equation by a power of moment and integrate them 
over all moment arguments. The zeroth moment equation is given by 

- j c(l, 2)d%d^p2 +0^1 ;;^P2ic(l, 2)d%d% = 0, (0th moment) (21) 
X = X21 = X2 — Xi, P21 = P2 ~ Pi- (22) 

As we can see, the time evolution of the zeroth moment depends on the first moment 
in the BBGKY equation. 

The first moment equation is given by 

- I c{l,2)pWpid'p2 + ^ I ^p-2lPUMd%d'p2 



+ / ^{Kl)c(2, 3) + d}'^d''x^d^p^d^pid^p2 
J a X31 

^{6(2)c(3, 1) + d}^d'xsd^p^d'pid^p2 = 0. 

Ci XqO 



^32 

(1st moment) (23) 
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This equation includes the second moment and then we need the second moment 
equation as follows: 

d I" ,^ ft ^ ,n ,5 d 



j c{l,2)pl,pl,d'p,£p2 + I ^p^,pl,p^2ic{l,2)d'p,d'p2 

^{6(l)c(2, 3) + d}{p\^- - p^^^--f}d'xsd%d%d% 

^ ^31 -^31 

- / ^{6(2)c(3, 1) + - Pi^}d'x,d'psd'p,d'p, = 0. 

J a X32 0:32 

(2nd moment) (24) 

As we can see above, the time evolution of the A'^-th moment depends on the + 1-th 
moment. So we should take assumptions in order to close these equations. DP used the 
assumption that the skewness of the velocity is equal to 0. In this paper, we do not assume 
about the skewness of the velocity fields in order to study the relation between the skewness 
and the two-point correlation function. 

Here we define the two-point correlation function ^, the mean relative peculiar velocity 
{v°') , the relative peculiar velocity dispersion H, E and the mean third moment {v°'v^v'^) as 
follows: 

n^a^^ = y c(2, l)d^pid^p2, (25) 
nV(l + Oma(^;") = J c{2,l)p^^d^pid^p2, (26) 
n'a%l + Oimaf[UP-^ + EPf] = J ci2,l)p^,,pWpid'p2, (27) 



r— T3 / P2P2lP2lP2ld^Pld^P2 

(ma) J 

P2d^Pid^P2 
(ma)%V(l + ' 



(28) 



where 

pf = pf = - (29) 



x"^ 



Here n is the mean density of the universe. The ^ which is defined in eq.(25) is equivalent 
to the two-point correlation function ^2{x) which is defined by another familiar definition 
as follows: 
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\ Mr) - {p)){p{r + x)- (p)))^^ . , ^ 



(P) 

= -^,Mr)p{r + x))-{p){p)\ 

- (P2(l, 2) - b{l)b{2))d'p,d'p2 

= ax2i). (30) 
where (p) is the mean density of the universe. 

In eq.(26), {v'^) is the mean relative pecuhar velocity given by 

/ P2{2,l)—p^,d'p,d'p2 

{v'^) = I — nio 

j P2{2, l)d^Pld^P2 

c(2, l)p^^,d'p,d% 



n^a^(ma)(l + 



(31) 



In eq.(27), 11 and E are the parallel and transverse correlated parts of the relative 
peculiar velocity dispersion of correlated particles, respectively. 

Here we define the skewness as follows: 



{{v-{v)r){{v-{v)r){{v-{v)r). (32) 



Prom the symmetry of the universe(homogeneity and isotropy), we can write 
(v") = {v)x°- /x. Hence, 



lyl^ ly"^ ( rt^^ ryl^ ry ] 

I •Aj •Aj / \ I I iQ I.C /,-., \ iX. 



(^;"^;V) = -2{vf \ + {v) I — (^;V) + —{v^v'') + —{v'^v^) \ + s"^^. (33) 



IX X X 



Furthermore the skewness can be written by the symmetry of the universe as follows: 



= n + «^ Ppu. (34) 



where the subscripts p and t represent the parallel and transverse component of the two 
particles, respectively. 
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Pppt and Put vanish because of the symmetry of the universe. 
Prom eqs.(28) and (33), we can get the following equations: 

J c(2, l)p^,pipl,d^pid'p2 = n-'a^l + 0(ma)3{[3(t;)(n - E) - 2{v) 

2{vj) 



(3 'y 



or. .^)(— '^^^ + — + — + ^'""^■y (36) 
3(1 + ^) XXX 



Then moment equations are following: 



dt 



(0th moment) 



(37) 



nV — [(l + Oma(T;0] + 



(l + 0[nP|f^ + EP°''] 



+ 



+ 



a 

Gm? 
a 

a 



-na 



/ c{2,2,)'^(fx^d'p^d'p2 



I c{3,l)^d^xsd%d^p, 

/ C(% - ^)d^x^ = 0, (1st moment) (38) 



9_ 
di 

+ 
+ 
+ 
+ 



n^a^(ma)^ 


d 




dx"" 


n^a^(ma)^ 


d 


ma? 


dx'^ 




d 



rftOt rf*l3 rf*'~^ 

3l 



(i + 0(^)(s + 



3(1 + a; X X ' 



-na 



[c{2,3M,^+p^2i^]d'x,d%d'p, 



+ 



Grrf 

a 

Gnf. 
a 



^31 

r.0 



-na 



-na 



' rc{3,iMi^+pi,^]d'x,d'p,d'p, 

^32 ^32 

' / rfbll(4 - 4) +P2l(# - ^)]d'xsd'psd'p,d'p, = 0. 
•> •^'31 -^32 -^31 -^32 

(2nd moment) 



(39) 
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where ( is the three-point correlation function defined by 

n^a\ = J dd?pid^p2d?p^. (40) 



2.3. Contraction of the Equation 

The zeroth moment and first moment equations are the time evolution equations of 
the ^ and (f), respectively. The second moment equation is the time evolution equation 
of the n, S. For convenience, we transform these equations by taking divergence of the 
first moment equation and by operating the following two operators to the second moment 
equation, 

— — = -(5^^ - (41) 

Hence we get two equations from the second moment equation. We call them 
contraction land contraction 2 equations hereafter, shown later in eqs.(50) and (51) , 
respectively. Here we assume the following relation, 

Cl23 = Q(ei263 + 63^31 + ^31^12). (42) 

Then the fifth term of the first moment equation ( |5B| ) is rewritten by 

OX'^ J X%i X%2 Ci UX'^ J X%i 



a 



a oxf^ J X31 

= 2^n'a'Q/^ [ ^^'^3(6263 + 63^31 + ^31^12) 

= 2^#a^Q A I ^rf^X3(62 + 61)63 
a ox''' 



where 



a ox'^ 

= 2 n^a^g— x + T 

a ox 

Grn? o 1 d , o^, 

= 2 n^a^g— — x^T 

a x^ ox 

= 2^n'a'QM\^[x'e], (43) 
a x^ ox 



X^T ^ X^^xfM ^ I ^{62 + e3l}63rf'x3. 



x = X2i. (44) 
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This definition is well defined for the vector component of x and x^i because of the 
symmetry of the universe. 



The seventh term in eq.(39) is rewritten by 

-I d[pl,{- 

'^-'31 -^32 -^31 "''32 

f X 



P P 7 7 

Ml / 31 -^32 \ , 3 /-^^l -^32x1 73 73 73 73 
/ d[P2l{-r - Zr) +^^2l(-r - -1^)]^ Xsd psd pid P2 



a dx/^dx'^ 

Gm^ 52 . 



= 4- 



.2 



^Gm', 



^(ma)n^a^g*^^^/|i[(^2\)ei263 + ((^la) + (^3\))636i + ^vl,)Ui^2\d^ x^ 
- 4^(ma)n a Q y ^[x ^^1263 + {(^ - ^ + ^ ^16361]^? X3 

= 4^(ma)nW^|^ /4^a;n^62 + ^e3i}e23ci^^3 
a oxf^ox^ J X31 X F ~ -^1 

3 ' 

1^31 

where 



+ 4^(ma)#a^g*^|^ / ^.^^ - ^^i.^i..d'x^ (45) 
a axPax^ J xi-, z \x — z\ 



^n^a^Cma{v2i) = f dp^^d^p (46) 

In general this relation (eq. [46]) is not satisfied. But DP showed the existence of d 
which satisfies this relation. Furthermore, in the strongly nonlinear regime, which we are 
interested in, we find from the dimensional analysis that this relation is correct in general. 

We can see from cqs.(39) and (45) that the seventh term of the second moment 
equation is transformed by operator A^''' as follows, 

^A^- / rfbli(4 - 4) - %d'xs£psd'pid'p2 



2^31 ^32 ^31 ^32 

^ 4^A^^ / dp^^i^d^'xsd'p.d^pid^p^ 

= 4^(ma)n3a^g*A/^- / #^^^62 + ^X.^]i.,d^x, 

+ 4^(ma)#a^g*A/^- / ^^^^ " ^16361^^X3. (47) 
a J X2,i z \x — z\ 

Finally we obtain the following four equations: 

■ d^ n^a^{ma) 1 d 
dt ma? x^ dx 



n^a'§- + ^'«'(^«) J_|_[^2(^ ^ ^^^^^1 ^ Q_ (Q^j^ moment) (48) 
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+ 



n^ar{maY 1 d d 2 



ma 



dx dx 



{ — (a;^(l + C)n)-2a;(l + C)S} 



a 



a x^ ox 



(1st moment) (49) 



d imaY d 



dt x'^ dx 
+ 



^[x2(l + 0n]-2a;(l + 0S 



+ 



+ 







n^a^{maY 3 


.9 1 




5x X 


n'^a^{maf 1 


d^ , 



x\i + omv){ii-i:)-2{vf} 



d 1 



+ 3(1 + 



[^'(l + e)«||] 



H na'^n^a^ — — b^fl + ^)(ma){v)] 

x^ dx 



a 

Gw? 



+ 4- 



a 

Gw? 



a 



{ma)rra^Q 
(ma)n^a^Q* 



dx^dx'^ J xli 
d' 

dxf^dx^ 



/ <-y»'-' rY* 



{V23) 



\x — z 

{V23) 



\x — z\ 



K23Ud''x3 = 



+ 
+ 



(2nd moment : contraction 1) 

n^a^imaY 1 d 



x'^ dx 



x 



rror{maY 119 
ma^ x"^ x^ dx 
Gm^ 



^(1 + 0(^){S 

.13 



3(1 + 



a J Xqi Xqi 



Gw/ 



^31 



^31 

7 



na'A^^ / c(3, l)b^i^ +pfil^]d^X3d3p3rf='pi 



a 

.Gm^ . .n q 

+ 4 {majn' a' ' 

a 

Gm'' 
+ 4 [ma]n 



- ™3 

•^32 



- ™3 - 
•^■32 



31 yJb A/\ 



{V23) 



(50) 



2^31 -yr(V3l) (^^23) 



X — ^ 



(2nd moment : contraction 2) 



}^23Ud'x3 = 0. 

(51) 
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We are interested in the strongly nonlinear regime and then the solutions of the above 
equations are expected to obey the power law because the self-gravity is scale-free. We 
assume that ^ is represented by the power law form as follows: 

e = ioa^x-^ (52) 

where l^i and 7 are constants. 

Then we obtain from the dimensional analysis in eq.(48) in the strongly nonlinear 
regime. 

{v) — —hax (53) 

where /i is a constant. 

Then the sixth term of the second moment equation (contraction 1, eq.50), is rewritten 

by 

A--{ma)n a Q J ^[-^^62 + ^CaiKd^ -3 

a oxPox ' J 



a 



maYrC'a'Q\-aK)^^{x^x^T\ 



= 4-=— -(ma)#a9g*(-d/i) — ^[,T^T] 
a x^ ox^ 

= A^{ma)n^a^Q*M{-ah)-^[x^e\- (54) 
The seventh term of the second moment equation(contraction 1, eq.50) is rewritten by 

4^(„„)*i3„,Q.^ 1 4,,{M _ <!^}&&,d3.3 = 0. (55) 

a oxPox^ J X31 z \x — z\ 

The sixth term of the second moment equation (contraction 2, eq.51) is rewritten by , 

.Gm^ ( n^3„9^*a/3t /■ ^31 r ^ r (^21) ^ , (^^23) 



4^(ma)#a^Q*A/^- / ^[^7{i!^ei2 + ^6iK23]d^X3. 
^ A^(ma)n'a'Q*(-ah)A^'^ / ^0:^62 + 61 K23rf'x3 

(X J Xqi 



.Gm^ 

''-^31 

a 



= 0. (56) 
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The seventh term of the second moment equation (contraction 2, eq.51) is rewritten by, 



a J xt 



X'^i z \x — z\ 



M23^3ld'xs = 0. 



(57) 



Finally we obtain the following four equations: 



(Oth moment) 



(58) 



1 9 1 5 r 2/1 /-N / \ 



+ STrGmna^ 

+ 2GmnaQM— — [x^C'^] = 0, 
x^ ox 

(1st moment) 



(59) 



1 9 2 1 5 



a? dt x^ dx 

1 1 93 



^[x2(i + e)n]-2x(i + os 



+ 



+ 



a X- 



dx^ 



x\l + 0{Hv){U-E)-2{vf} 



ax"^ dx X dx 



ax"^ dx X dx 



4 a 1 



1 (9 

+ 16nGmna—^—x^(l + E)(v) 
x^ dx 

1 52 

+ 4GmnaQ*(dh) — — [x'^M^^] = 
x^ dx'' 

(2nd moment : contraction 1) 



0, 



(60) 



1 d 



1 1 d 



a^dt 



+ 



a x^ dx 
1 1 d 



a x^ dx 
+ J = 



x\l + ^){v){T.+ 



3(1 + 



(2nd moment : contraction 2) 



(61) 
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where, 

J ^ ^na'A^^ I c(2, 3)^^ + p'2i''^]d'x,d'p,d'p, 
a J x^i x^i 

- ^na'A^^ [ c(3, IM,^ + pi^d'x.d^d'p^ (62) 

a J X^2 ^32 

This term is negligible in the strongly nonlinear regime. In deriving the above equations, 
we made use of the symmetry of the background universe and assumed that the three-point 
correlation function can be written by the products of the two-point correlation functions. 
And also ^ is assumed to be given by eq.(52) which is expected to be correct in the strongly 
nonlinear region. 



3. SCALE-INVARIANT SOLUTIONS IN THE STRONGLY NON-LINEAR 

REGIME 

In this section, we discuss how the BBGKY equations are approximated in the strongly 
nonlinear limit. And we discuss how the power index of the two-point correlation function 
relates with the skewness, three-point correlation function and the mean relative peculiar 
velocity. Furthermore we investigate whether the stable condition is correct or not. 



3.1. BBGKY Equations in the Nonlinear Limit 

In the strongly nonlinear regime, the two-point correlation function is much larger than 
unity, ^ 1 at X <^ 1. In this limit, the BBGKY equations are given by 



at a x^ ox 



(0th moment) 



(63) 



1 1 d , d 



1 d 



a x"^ dx dx 



{ — (x^^U) - 2x^J:} + 2GmnaQM— — \x^^^] = 0, (1st moment) 



(64) 



a? dt x^ dx 



d_ 
dx 



[x^^U] - 2xiE 



1 1 



13 8 18 
ax"^ dx X dx 



.i^f {3{i;)(n - E) + - 3sx)} 



-17- 



1 (9^ 

- iGmnaQ*(dh) — — \x^M^^] = 0, 



(2nd moment : contraction 1) (65) 



0. 



(2nd moment : contraction 2) (66) 

These equations are equivalent to eqs.(39) ^ (42) in Ruamsuwan & Fry (1992) while 
some notations are different from each other. The zeroth moment equation (63) is derived 
without any assumption. Here it must be noted, however, that there is assumption about 
the three-point correlation function in deriving the first moment equation (64). And in the 
second moment equations (65) and (66), there are the terms of the skewness and so we 
usually need the higher moment equations in order to solve eqs.(65) and (66) while we do 
not need them in our analysis. 



3.2. Scale-Invciriant Solutions in the Strongly Nonlinecir Limit 

In the strongly nonlinear regime, it is naturally expected that the effect of the 
nonlinear gravitational clustering dominates and then the solutions in this regime have no 
characteristic scales, that is, they arc expected to obey the power law due to the scale-free 
of the gravity. So wc investigate the power law solutions of the ^, {v),Il and E. Then we 
assume that the two-point correlation function ^ is given by 

e = ioaf'x-'. (67) 

We can see from eq.(63) that the mean relative peculiar velocity (v) is given by the 
dimensional analysis as follows: 

(v) = -hdx (68) 

/5 = (3 - ^,')h (69) 

The stable condition means that h = 1 because (f) = ax + (ax) = ax + (v) = 0. 
Hereafter we call this parameter h, relative velocity parameter. We obtain the power law 
solutions of the other quantities from eqs.(64) and (65) as follows: 

H = Uoa^-^x^-^, (70) 

E = Eoa'^-'a;'-^ (71) 

S|| = siioda'^-^x^-^, (72) 

Si = si_oda^-^x^-^. (73) 
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Here Ho, Sq, syo and s^o ^-re constants. 

Prom the above results, we obtain the next relation from eq.(66). 



(1 + 2/3) -(7 



27)(/i- A) = 0. 



(74) 



Here A is defined by 



A 



(75) 



So 



We can find from eqs.(69) and (74) that 



= 1 + (7 - 27) A. 



(76) 



Here it must be noted that eq.(76) is determined on the assumption that the three-point 
correlation function can be written by the products of the two-point correlation functions(see 
eq.[42]). Then if this assumption is correct and also the skewness vanishes, which means 
A = 0, we can see from eq.(76) that the relative velocity parameter h is equal to unity. This 
fact means that the stable condition h — 1 is not an assumption, but should be satisfied in 
the strongly nonlinear regime when the three-point correlation function can be represented 
by the products of the two-point correlation functions and the skewness vanishes. However 
it may be physically natural that the skewness does not vanish in the nonlinear regime 
while the skewness equals in the linear regime. Then h should not be equal to 1. This 
means that the stable condition is satisfied for the novel case that the skewness vanishes. 
Furthermore h varies when the assumption about the three-point correlation function ( 
changes. For example, we assume that C c< C^^^^^^ here 5 is a constant. Then we can find 
from eqs.(63) ~ (65) that 




(77) 
(78) 
(79) 
(80) 
(81) 
(82) 



Prom eq.(66), it is found that 



2/3(1 + S) + l-{7- 27(1 + S)}{h - A) = 0, 



(83) 



and then 



l + {7-27(l + (5)}A 
1-65 



(84) 
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Thus we can see that the value of the relative velocity parameter h also depends on the 
assumption about the three-point correlation function. It must be noted here that a little 
change of 6 results in the large change of h due to the factor (1 — 6(5) in eq.(84). 



In the previous subsection, we found that the relative velocity parameter h depends on 
the three-point correlation function and the skewness. Then, in this subsection, we discuss 
about the allowed value of h from the physical point of view. 

In general, the mean relative velocity (r) is represented as follows: 



DP assumed the stable condition in which h = 1 in the virializcd objects. But wc found 
that h depends on the ^ and A. The stable condition is satisfied only in the novel case that 
A = if eq.(42) is correct. We show here that even in the virialized region h can have other 
values rather than unity. For example, we consider the spherical collapse for simplicity. In 
this case, it is found that a object has collapsed and virialized when the linearly estimated 
density fluctuation Sun of this object becomes just 1.69(= 5c) (Peebles 1993). When the Sun 
of the virialized object with mass scale M becomes greater than Sc, then the object with 
mass scale M' > M should collapse and virialized at the same place(Press & Schechter 
1974; Yano, Ngashima & Gouda 1996). This means physically that matters around the 
virialized object accrete to this object with mass M and the larger virialized object with 
mass M' are newly formed. In this situation, the linearly estimated density fiuctuations 
of the virialized object is given by Sun = Sc, which means the real density fluctuations 
Pnoniin = IStt^ (p) (Suto, 1993). Here (p) is the mean density of the universe at the collapsed 
time. Thus the density of the virialized object Pnoniin is in proportion to (p) oc a~^, where 
Qc is the scale factor when the object has just virialized. 

Then we consider a sphere of the virialized object with physical radius r and consider 
the change of the mass of the shell at r. The mass change is represented as follows: 



3.3. Relative Velocity Parameter 



(r) = {ax + ax), 

and then the mean relative peculiar velocity (v) is given by 

(v) = (ax) 

— —ax + (r) 
= —hax. 



(86) 



(85) 
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4 o, , d 
-.r'(-3)-p. 



(87) 



Together with r = ax, we obtain 



(f) 



= ax. 



8. 



.8 



And then the mean relative pecuhar velocity is given by 

(f ) = —ax + (f) 



= 



(89) 



In this case, the relative velocity parameter h is equal to 0. This case is one example 
of the physical situation that h is not equal to 1 even in the virialized regions. 

In the hierarchical clustering picture, the variance of density fluctuations on small 
scales is larger than that on the larger mass scale and the smaller objects are formed at 
first and after that the larger objects are formed. In this case, we can consider two extreme 
physical situations in the process that the larger objects are formed. One of them is as 
follows; the smaller collapsed objects have clustered and survived in the larger object. In 
this situation, the objects which are formed are never broken and the mean interval of 
particles in the objects does not change. Hence the relative velocity parameter h is equal 
to 1. Another situation is as follows; matters around the virialized object accrete to this 
object and the large virialized object are newly formed. In this situation, the mean interval 
of particles expand with the Hubble velocity as shown in eq.(88). 

From the above argument, we can conclude that the way of clustering decides the value 
of the relative velocity parameter h. We believe that the real way of the clustering might 
be between the above two extreme situations. The relative velocity paraemter h should not 
be greater than 1 because (r) have to be positive. And as we mentioned above, h — is 
the lower limit of h for the virialized object. Thus the relative parameter h takes the value 
between and 1, 



This condition gives constraint on the skewness A and the three-point correlation function 



0<h<l. 



(90) 



C- 



3.4. 



Self-Similar Solutions 



DP showed the existence of the self-similar solutions under the stable condition and the 
assumptions about the three-point correlation function and the skewness. In the self-similar 
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solutions, the power index of the two-point correlation function in the strongly nonlinear 
regime is related to the power index n of the initial power spectrum, where n is defined by 

P{k) (X A;". (91) 

In the hnear regime, the two-point correlation function is given by(Peebles 1980,1993), 

e oc a2x-(3+") oc (— )-(3+-) (92) 

where ^ 

a = . (93) 

3 + n ^ ^ 

On the other hand, in the strongly nonlinear regime, the two-point correlation function 
obeys the following evolution equation(Oth moment, eq.[63]). 

Here we used the relation (v) — —hax. 

This equation is rewritten by transforming the variables, x and a to the scaling variable, 
s = x/a^, 

ds~ P + hs 

If the self-similarity is satisfied for the scahng variable s, /3 should be equal to 
a = 2/(3 + n). 

Then it is found from eq.(95) that 

^ oc s"^ = s~^, (96) 

and the power index 7 is given by 

-,_ _ 3M3 + n) 

^ a + h 2 + h{n + 3)' ^ ' 

If the relative velocity parameter h — 1 is satisfied, 7 = 3(3 -|- n)/(5 -|- n), which is the 
result shown by DP. 

In the previous subsection, the allowed range of h is shown (0 < h < 1). Then we get 
the range of the power index 7 as follows: 

3(3 + n) , , 

- ' - 2 + n + 3 ^ ^ 



-22- 



Now we investigate the possibility that there exist the self-similar solutions in which 
the power index 7 is independent of the initial power index n. Wc can see from eq.(97) that 
this situation occurs when the relative velocity parameter h satisfies the following relation 
( Padmanabhan (1995)), 



where c is a constant. 

However h should be within the range between and 1. If c is not equal to 0, /i becomes 
the value of larger than unity for some n even if h is within the allowed region for other n. 
Hence only when c is equal to 0, that is, /i = 0, we can get the self-similar solution in which 
the power index 7 is irrespective of n and also h is always within the allowed region. This 
result is reasonable from the physical point of view shown in §3.3. When /i = 0, the smaller 
collapsed objects become to be absorbed hierarchically in the larger object. In this case, 
the information of the initial memories are erased. Then the power index 7 is irrespective 
of n. On the other hand, when h ^ 0, substructures in the larger cluster survived and the 
initial memories are not completely erased. 



4. RESULTS AND DISCUSSION 

In this paper, we investigated the scale-invariant solutions of the cosmological BBGKY 
equations in the strongly nonlinear regime. The mean relative peculiar velocity depends on 
the skewness and the three-point correlation function. The stable condition which DP used 
is satisfied for the novel case that the skewness vanishes and the three-point correlation 
function can be represented by the products of the two-point correlation functions. 
Furthermore the power index of the two-point correlation function in the strongly nonlinear 
regime depends on h, the skewness and the three-point correlation function(see eqs.[76] and 
[84]). 

In the hierarchical clustering, there are two extreme situations in the way of the 
clustering. One of them is that the collapsed object can not be broken and clustered 
together to form the larger cluster. In this situation, the mean separation of the particles 

does not change as time increases and h = 1. This corresponds to the stable condition. 
Another situation is that matters around the collapsed object accrete to this object and the 
larger collapsed object is newly formed. In this case, the mean separation of the particles 
is expanding with the Hubble velocity and then h = 0. It is believed that in general h 
should take the value between and 1. This constrains the skewness and the three-point 
correlation function. In general, the skewness does not vanish in the nonlinear regime and 
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so it might be impossible for the stable condition(/i = 1) to be satisfied. 

The self-similar solution is shown by DP under some assumptions and the stable 
condition. If there exist the self-similar solutions under the other conditions, then the power 
index 7 is given by 

0<7<^^, (100) 

because < h < 1. 

In general we may expect that the power index 7 does not depend on n because the 
systems forget the initial memories on the nonlinear regions due to the nonlinearity of the 
gravity. If h equals c/ (n + 3), then 7 is irrespective of n and is given by 3c/ (2 + c). However 
h should be between and 1. Then c should be zero, that is, h = 0. In this 
mentioned before, the collapsed objects are hierarchically absorbed into the larger objects 
and the substructures in the cluster can not survive. Then the initial memories are erased. 
Thus it is physically reasonable that 7 does not depend on the initial condition n in this 
case. 

In general, the systems have substructures in some clusters and no sub clusters in other 
clusters and so 7 depends on n. However the relation between 7 and n may not be given by 
the result shown by DP because the skewness does not vanish in general on the nonlinear 
regions and so h should not be equal to 1. Here we comment that 7 = is expected from 
the catastrophe theory. And we can see that if the solutions of the correlation function 
have the self-similarity, the power index 7 = 2, which Saslaw(1980) predicted, is impossible 
to appear stably in the strongly nonlinear regime irrespectively of the initial power index n 
because in this case c should be eqaul to 4 and. then /i > 1 for n < 1. 

Ruamsuwan & Fry tested the stability of the scale-invariant solutions derived by 
DP and found that they are marginally stable. We will analyze the instability of the 
scale-invariant solutions shown in this paper under the other assumptions beside those 
taken in DP. Furthermore we will investigate whether there exist the self-similar solutions 
under the various assumptions about the three-point correlation function and the skewness 
in the future. 
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